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Abstract 

We present a new type of modified Galerkin method. It is a construc- 
tion with several (inductively defined) levels, that provides approximate 
solutions of increasing accuracy with every new level. These solutions are 
constructed as approximations of the so called induced trajectories (no- 
tion on which the definition of a class of approximate inertial manifolds 
used in the nonlinear and postprocessed Galerkin methods is based). 
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1 Introduction 

We consider the Navier-Stokes equations for a planar flow, with periodic bound- 
ary conditions. The functional framework and the basic hypothesis are presented 
in Section 2. 

As usual in the Galerkin method, the space of functions is split into a di- 
rect sum of two subspaces: one is the finite dimensional space spanned by the 
eigenfunctions corresponding to a finite set, T^, of eigenvalues of the linear 
operator A = — A and the other is the orthogonal complement of the first. 
The solution u of the Navier-Stokes equations will be projected on these spaces: 
u = Pu + Qu = p + q, where P is the projector on the finite dimensional space, 
and Q = I - P (Section 3). 

In Section 4 we improve the estimates for q proved in 3 . There, the 
[L^ (fj)] norm of this function is found to be less than KqL^S, with S = j^, 
where A is the least eigenvalue of A, A is the least eigenvalue of A not belong- 
ing to Fm, and L depends increasingly on the number of eigenvalues in r„i. The 
presence of L is not convenient for our work, since we construct an iterative 
approximation processus, and if we would use this estimate, at every step a 
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factor of i^/^ would appear in the evaluation of the error. This would lead 
us to bad estimates of the accuracy of our approximate solutions. We obtain 
estimates of q independent of L. 

Our modified Galerkin method is presented in Section 5. The first level of 
the method is related to the already classical postprocessed method Hj. This 
one consists in correcting the Galerkin approximation of the solution, computed 
at the end of the time integration interval, po (T) (that is, the solution of (|24p 
at T), by adding to it the function qo (T) = <&o(Po {T)). $o(0 is the function 
whose graph is the approximate inertial manifold (a.i.m.) defined in [3]. The 
approximation of u (T) is taken as po (T) -\-c\q{T). Unlike this "postprocessed" 
Galerkin method, in our method the function qo {t) — <I>o(Po (0) is computed 
at every moment t. From the numerical point of view this means that it must be 
computed at every point of the time grid on [0, T]. The approximate solution at 
every t is \io{t) — po {t) +qo(t). The error (in [Lpg^ (^)] ) of this approximation 
is of the order of 6^^^. We must remind here the notion of induced trajectory 
of [15j . There a family of functions is defined, {u^ j„}^.^q, that approximate the 
exact solution of the Navier-Stokes equations. The first of these is Uo,m(i) = 
p(t)+qo,m(Oi where qo,m(0 = 'I'o(p (0): ^-^id p(<) is, as above, the P projection 
of the exact solution. So, our function Uo(t) is an approximation of Uo,m(t). 

At the second level we look for a new (and better) approximation pi of p, by 
solving an equation closer to the P projection of the Navier-Stokes equation 
than the Galerkin equation. That is, in the equation for pi, in the argument of 
the nonlinear term, we approximate p + q (that appears in the projected N-S 
equation p3)) ) with Pi+qg, with qo defined above (equation ([27|) in Section 5). 
The initial condition for pi is Puq. The equation (P7)l is different from those 
arising in the non-linear Galerkin methods, since here qo is already determined 
at the preceding step, while in the non- linear Galerkin methods t2j in the 
equation for p, the argument of the nonhnear term is p-|-$o(p)- Solving equation 
(|27|) of our method is not essentially more difficult than solving the Galerkin 
equation. ^ ^ 

Then we compute qi (t) = <l'i(pi {t) , qo (t)) where <&i is given by ([28|) . The 
definition of qi (.) is inspired from that of the function qi.m (.) of [15] and qi (.) is 
an approximation of this latter function. There is an obvious connection between 
the definition of qi {t) and the second a.i.m. from the family defined in [15] and 
used in the non-linear Galerkin methods (the definition of this a.i.m. uses that 
of qi,m (•))■ The new approximate solution of the Navier-Stokes equations we 
define is Ui [t) = p^ (t) (t) . This is a better approximation of the exact 
solution that the preceding one, since the error is of the order of S'^' (as IS 
shown in Section 6). 

We define inductively the next levels of our method. By assuming that 
we already found pj, qj for j < fc + 1 (fc > 0), we construct the equation 
for Pk+2, taking Pk+2+(lk+i the argument of the nonlinear term (equation 
([30])). The "small eddies" component of the solution will be approximated by 
qfc+2 {t) = $fe+2(pfc+2 (t) ,q_k+i (t) ,q_k (t)), where $fc+2 is a function (given by 
the right hand side of (|3T|) ) whose construction was inspired by that of the 
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function qfe+2,m {t) of [H]. The error of Uk+2 (t) — Pk+2 (0 +^k+2 (0 is of the 
order of (5^/'*+'^'^+^)/^. The set {uk+2 {t) ', t > 0} is an approximation of the 
induced trajectory {ufc+2,m (t) ; ^ > 0} . 

In Section 6 we prove the estimates of the error of our approximate solutions. 
Finally some comments on the advantages and drawbacks of this method are 
given in Section 7. 

2 The equations, the functional framework 

The plane flow of an incompressible Newtonian fluid is modelled by the Navier- 
Stokes equations: 

f, (1) 

0, (2) 
uo (•) , (3) 

where u — u{t,x) is the fluid velocity, x efi C M'^, u (., x) : [0, cx)) R^, 
p (., x) : [0, oo) ^ M is the pressure of the fluid, v is the kinematic viscosity, and 
f is the volume force. We take here = (0,/) x (0,1) and consider the case of 
periodic boundary conditions. 

The kinematic viscosity is measured in centistokes (= mm? /s). For water 
at 20 Celsius degrees it's value is around 1. In order to have coherent measure 
units, we consider the velocity measured in mm/s. We do not focus here on the 
mechanics of fluids aspects of the problem, but we focus on the mathematical 
construction of the approximate solution. However, we must remark that the 
method is appropriate for the study of Newtonian fluids not having very small 
kinematic viscosity. 

We assume that f is independent of time and is an element of [L"^^^ (5^)] . 
As is usual in the study of the Navier-Stokes equations with periodic boundary 
conditions, we assume that [13], [12] 

f =i|^f(x)dx = 0, (4) 

and that the pressure is a periodic function on 17. For simplicity we will assume 
also that the average of the velocity over the periodicity cell is zero. 

The velocity u is thus looked for in the space 
H = |v; V e [Lpg^ (il)]^ , divv = 0, u = o|. The scalar product in Ti. is 

(u, v) = {uivi +U2V2)dx, (where u — (1*1,1*2) , v — (wi,W2)). The induced 
norm is denoted by | |. 

We also need the space V = |u e [Hp^^ (fi)] ^ , div u = 0,u = o| , with the 

scalar product ((u, v)) = j=i (ff^' ' ^'^"^ induced norms, denoted by 
II II . We denote A = -A and observe that A is defined on D{A) = HnH^ (ft) . 



du 
'dt 



uAu + (u • V) u+Vp = 

divu — 
u(0,.) = 
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We shall focus on finding approximations for the function u. 
The classical variational formulation of the Navier-Stokes equations [T3l leads 
to the abstract equation 

^-z/Au+(u- V)u f inV', (5) 
u(0) = uo, UoGW. (6) 

The notations B(u, v) = (u • V)v, B(u) ~ B(u, u), b(u, v, w) = (B(u, v), w) 
will be used below. 

For the bilinear application B(u, v) the following inequalities 

|B(u,v)| < ciluplAupllvll, (V) u ei?(A), V eV, (7) 

, (V) u eD(A), V eV. (8) 

hold [B], [in], [H]- We remind the following properties of the trilinear form 
b(u, v,w) (valid for periodic boundary conditions [12]): 

b(u,v,w) = -b(u,w,v), (9) 
b(u,v,v) = 0, (10) 

as well as the following inequalities [12] 

|b(u,v,w)| < C3|un|up||v|||wn|wp, (V)u,v,weV, (11) 
|b(u,v,w)| < C4|u|^||u||^||v||5|Av|^|w|, (V)ueV, veD(A),we7^. (12) 

For the problem ([5]), ([6]) we have the classical existence and uniqueness 
results for the equations Navier-Stokes in R^, with periodic boundary conditions. 

Theorem 1 [13]. a) Ij uq G H, f G 7i, then the problem ([5]), ([6]) has 
an unique solution u G C'' {[0,T];TI) n (0, T; V) . b) If, in addition to the 
hypotheses in a), Uo G V, then u G {[0,T];V) (1 {0,T; D{A)) . The 
solution is, in this latter case, analytic in time on the positive real axis. 

The semi-dynamical system {S (^)}t>o generated by problem ^ is dissipa- 
tive [2], [H]. More precisely, there is a Po > such that for every i? > 0, 
there is a to{R) > with the property that for every Uq G 7i with |uo| < R, 
we have |S'(t)uo| < Pq for t > to{R) . In addition, there are absorbing balls 
in V and D (A) for {S{t)}^yQ, i.e. there are p;^^ > 0, p2 > and ti{R), t2{R) 
with t2{R) > ti{R) > to{R) such that for every R > 0, |uo| < R implies 
115 {t) uoll < Pi for t > ti{R) and |A5 {t) uo| < for t > t2{R). 



|B(u,v)| < C2||U||||V|| 



1 + ln 



iAur 

A, Hull 
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3 The decomposition of the space, the projec- 
tions of the equations 

The eigenvalues of A are A^.j, = ^ (j? + j|) , (ji, ^2) e {(0,0)} , and 
the corresponding eigenfunctions are 



I \i\ V I 



"^juj. I |j| --y- I 

1 

where |j| = (^jf + j|) ^ [15]. These eigenfunctions form a total system for Ti.. In 
order to be able to write easily sums involving the four eigenfunctions above, 
we denote them as follows 

s+ 1 s— 2 c+ 3 c— 4 

For a fixed m e N we consider the set r,„ of eigenvalues Ajj.j^ having 
< ji, j2 < TO- We define 

A ■ -A -A 

47r2 2 
A : = A,„+i,o = Ao,TO+i = (m + 1) , 



(5 = (5(m):=^ = 



2 ■ 



A (m+l) 

A is the least eigenvalue not belonging to r„j. The eigenfunctions corre- 
sponding to the eigenvalues of span a finite-dimensional subspace of Ti. We 
denote by P the orthogonal projection operator on this subspace and by Q the 
orthogonal projection operator on the complementary subspace. We write for 
the solution u of ([5]), 

p = Pu, q = Qu. 

By projecting equation ([5]) on the above constructed spaces, we obtain 

^-i.Ap-^PB(p-fq) = Pf, (13) 
^-i.Aq+QB(p + q) = Qf. (14) 
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4 New estimates for the "small" component of 
the solution 

In [3] is proved that for every R > 0, there is a moment ts (i?) > t2{R) such 
that for every |uo| < R, 

|q(<)| < KoLU, \\q{t)\\ < K^lUK (15) 
|q'(i)| < KlU, |Aq(i)|<X2i^ t>t3{R), 

where Kq, K'q, Ki, K2 depend of i^, |/| , A and, for the way we chose the 
projection subspaces, L = L{m) = 1 + ln2m^ (see also [U]). The constant L 
comes from the use of inequahty ([5]) in the course of the proof of (|15p . More 
specific 

L = sup I 1 + In — — i^-'-TT I = max I 1 + In — 




Aer,„ V Ai 



In2m2 



In the sequel we shall improve the above estimates, trying to eliminate L 
(which tends to infinity with m) from the constants. The idea is that of refining 
the contribution of the term QB(p) resulting from QB(p + q) in p4|) . We start 
from the trigonometric relation 



. , jlXi±j2X2\ . kiXi±k2X2\ 1 

sm I 2tt : I sm | 27r -, | = - 



o (Ji ~ ^1)2^1 ± (j2 - fc2)a;2 

cos ZTT J 

(jl + ki)xi± {j2 + k2)x2 



cos 27r 



I 



and the similar ones for all other combinations of sine and cosine that might ap- 

4 

pear in the scalar product of two eigenfunctions. Since p = ^ ^ p*^ j2^ji j2 ' 

0<iij2<ni»=l 

from(pV)p=( J2 EPn.n'^n,],^] { E E PL,fc2Wfei.fe2 ) ' only 

\0<jij2<m«=l / \0<ki,k2<ml = l j 

those products of terms that have ji + fci > m + 1 or j2 + fc2 > rn + 1 will belong 
to QW. 

We consider from this point on that to is even, and we set m = In. 

If for wj^ and wj,^ we have ji, j2 < n and fci, fc2 < then 
(^ii 52^) ^fci fc2 belongs to P7i. We are led to the idea of considering the 
subspace of 7i spanned by all the eigenfunctions w*^ with < ji, j2 < 
1 < i < 4. We denote by Pp the projection operator on this space and set 
Pq = P — Pp, Pp = Ppp and pg = Pqp. Obviously 

Q (PpV) pp= 0. 
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|q(t)l 


< 


Co5, 


l|q(t)l! 


< 


CiS^ 


|q'(t)l 


< 


CqS, 


|Aq(t)| 


< 


C2, 



On another hand, we see that is a truncation of (I — Pp)u, hence 
|p,| < 1(1 - Pp)u| , ||p,|| < 11(1 - Pp)u|| . Then, by setting Si ^ S (n) = 
Li = L{n) = 1 + ln2n^, the estimates (IT5|) imply 

|p,| < KqL15i, llpgll <ifiif4: (16) 
|p;| < K'^LjSi, |Ap,|<i^24. (17) 

We use these inequahties in order to refine the estimates (|T5|l . In the rest of 
the paper we shall assume that for a fixed i? > 0, the function Uq is such that 
|uo| < R. We state and prove 

Theorem 1. There are some constants Cq, Ci, Cq, C2, depending only 
on V, A, |Qf| such that, for t large enough, the inequalities 

(18) 
(19) 
(20) 
(21) 

hold. 

Proof. We have, with the notation settled before the Proposition, 

QB(p) = QB(pp+p,) = QB(pp) + QB(pp,p,) + QB(p„Pp) + QB(p,) 
= QB(pp,p,) + QB(p„Pp) + QB(p,), (22) 

since QB(pp) — 0. Now, as is usual, we take the scalar product of (|14p with q, 
and by using (fTU]) and (|22p . we obtain 

^^ + ^^'1'^"' - l(B(Pp,P,),q)| + |(B(p„p,),q)| + 

+ |(B(p,),q)| + |(B(q,p),q)| + |(Qf,q)| . (23) 

For the first term of the right-hand side, the following estimates (obtained by 
using and (Ull)) hold 



|(B(pp,p,),q)| < Ci|p,|^/^|Ap,|i/2|lp,|l|q|<CipyV^/^ifi44-i^ ||q|l 

< c?PoP2i^fii'5lJ^ + ^l|q||'- 
For the second term we obtain the inequalities 

|(B(p„Pp),q)| < ci|p,|i/^|Ap,|i/^||p,|||q|<cii^y%^(5fp2^'pi-^ i|q|| 

< clK,plp,L\'h,^ + '^\H\\ 
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For the third term we have 

|(B(p,),q)| < C2Li\\p,\f\q\<C2L-^KfLi6i^\\q\\ 

A2 

< clKtLLlSl^ + !^M\ 
and for the fourth, by using (I12p 

|(B(q,p),q)| < C4|qp||qp llpplApHql < 

< cK'.^'Lis'^Kl^'Lh'^plpl^ ||q|| 

< ^clKoK.p.p^LS^ + ^\\cif. 

At last ^ 

|(Qf,q)|<|Qf||q|<^^ + ^||qf . 

The above inequahties and (|^ lead us to 

^ d 2 , ol'A 2 , ^2r 

with 



- — 



clpoP^KfLiSi + cIKipIp^lI'^Si + clKtLLl5l + 
+clKQKip^P2LSi + IQff 
It follows, with the usual Gronwall Lemma, 



|q(t)|^<|q(0)re-i''^* 



8C0 r2 



f i/At 



< i^'^' for 



hence, for i4 (i?) > t3{R), taken as to have |q(0)|' 
t > ti (R) , we obtain with Co = 

The functions of n : LiSi ^ L (n) 6 (n) (l + In 2n'^) / {n + if , 
lY^Si = ^/L{n)5 {n) = ^/TThPh? / {n + if, 

LLlSl = L{2n)L{nf 5{nf = (l + lnSn^) {l + \Ti2n^f / {n ^ if and 

L5^ = (1 +ln8n^) / (2n + 1)'^, that appear in the structure of Cq, have at 
n = 2 values less than 1 and are decreasing when n increases (for n > 2). Then, 
for n > 2 

C^o < ^ {cIpoP2K! + cIK^pIp2 + clKf + cjKoK.p.p^ + |Qf |') 
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and the right hand side depends only on A, |Qf | . More than that, since all 
the functions defined above tend to zero when n — > oo, we can choose n large 
enough so that ^^^^^ becomes the dominant term in Cq. 

For that n, Co will be of the order of . However, the structure of 
Kq, Ki, Pq, P2 show that if ly is very small, then n with the above 
property must be very large. 

Now, we aim to estimate ||q|| . By multiplying equation by Aq, and by 
using we obtain 

^^^^ + '^|Aq|' < |(B(pp,p^ + q),Aq)| + |(B(p,+q,Pp),Aq)| + 

+ |(B(p,+q,Pg+q),Aq)| + 
+ l(Qf,Aq)|. 

For the first term in the right hand side we have ([7]) 

|(B(pp,p, + q),Aq)| < Ci|ppp|App|^||p, + q|||Aq| 

< cpI^'pI^'^K.sI |Aq| 

< clpoP2KfLi5i-+'^\Aq\\ 



V 



for the second, with ([T2 



|(B(p,+q,Pp),Aq)| < C4|Pg+q|5 ||p,+q||3 ||pp||2 |App|2 |Aq| 
< c^ClblLlKlblplpl |Aq| 

V o 

for the third, also with (|12p . 

|(B(p,+q,Pg+q),Aq)| < C4 |pg+qp||p,+q|l | A (p,+q) | ^ | Aq| 

< C4c|(5f LfXifJfL^-ft:! |Aq| 

V o 

and for the fourth we have 

|(Qf,Aq)|<^ + ^|Aqr 

V 8 



We denote 



V V 
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Then the differential inequality for |lq|| becomes 



d ||q|| 

dt 

that yields 



^A||qf <C2, 



||q(t)ir<||q(0)|re-^^* + i-C?<S. 

Let t4 (R) such that for t > [R) > (R) the inequality |lq (0)||^ e"'''^* < 
^CfS holds. For t > (R) ^ holds with Ci = ^J^Cl. 

We remark that Li5i ^ L (n) 6 (n) , Lf6J = L (ny^ S (n)^ , LjSl = 

- 2 

L{n)^ S (n) have values less than 1 for n — 2, decrease when n increases for 
n > 2 and tend to zero when n — > oo. Hence, Ci may be replaced with a 
coefficient that depends only on i^, A, |Qf | and not on n. 

Moreover, for n large enough, each of the first four terms of Ci becomes 
smaller than ^^7=. 

As for the solution u in [T3] , it can be proved that q (t) is analytic in time 
and is the restriction to the real axis of an analytic function of complex variable 
defined on a neighborhood of the real axis, and by using the Cauchy formula, 
we obtain (|20|) . 

Finally, from (fT4|) we have 



V 



^ + QB(p + q)-Qf 



and with the above estimates we obtain (1211). □ 



5 The new modified Galerkin method 

Let us fix a T > ^4 (i?). The interval [0,r] is the interval on which we seek 
the approximate solution. Obviously, all the above inequalities are valid for 

t e \u (R) , T] . 

In this section we just present the method, while in the following section we 
estimate the error of the method. 

5.1 The first level 

The first level of our method is related to the post-processed Galerkin method 

of a. 

Let po {t, x) be the solution of the equation (the Galerkin approximation of 

my- 

p^, - :.Apo + PB (po) = Pf, (24) 
Po(0) = Puo, 
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and 

qo(t)= *o (poW) , 

where $o ■ QT^ is the function whose graph is the a.i.m. Mq defined in 

[3], that is 

*o(p)-(^^A)-MQf-QB(p)]. (25) 

We define the corresponding approximate solution for the Navier-Stokes 
problem ([S])-® as 

uo (t) = po (t) + qo (t) . (26) 

Unhke the method of [3], we compute qo at every moment of time, and 
not only at the end of the time interval, T (in the course of the numerical 
implementation of this method, qo will be computed at every point of the grid 
on [0,T], constructed for the integration of ([24]) ) . 

Remark. In [15] the function qo,m {t) = *o (p (t)) is defined (with p {t) the 
P projection of the exact solution), and then the function Uo,m {t) = p{t) + 
qo,m (t) , is constructed, it's positive trajectory being named "an induced tra- 
jectory". Since po (<) is an approximation of p (t), as is proved in Section 7, it 
follows that {uo {t); t > 0} , is an approximation of this first induced trajectory 
of [H]. 

5.2 The second level 

The next level is different from both the nonlinear Galerkin methods and the 
post-processed Galerkin method already defined in literature. At this level we 
make use of qo calculated in the previous step and define pi as the solution of 
the equation: 

p'i-j.Api+PB(pi+qo) = Pf, (27) 
pi(0) = Puo. 

We expect this correction of po to be closer to p than po itself. Then we set 

qi(t) = (i/A)-MQf-QB(pi(i))-QB(pi(t),qo(i)) - 

-QB(qo(i),Pi(i))]. (28) 

We define the approximate solution for (I5|)-(l6|) at this level by 

ui(t) = pi(i)+qi(i). (29) 

Remarks. 

1. In the non-linear Galerkin method for the approximation pi of p, an 
equation, similar to (I27p . but with PB (pi+4>o (pi)) instead of PB (pi+qo) is 
considered. 

2. In what concerns qi, the right hand side of ([28]) (let us denote it by 
$i(pi (t) , qo (t)) ) was inspired from the function qi_.m (t) of [15] . 
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This one is defined as 

qi,™(t) = (z.A)-i [Qf-QB(p(<))-QB(p(t),qo,„(<)) - 
-QB(qo,™(0,pW)] 

and with it's help the function Ui^™ (t) = p{t) + qi.m (Oj defined, that 
generates a new induced trajectory. The construction of the second a.i.m. in 
[15], Ml, is based upon the definition of function qi.m (t) . Mi is the graph of 
a function $i : PH QH, given by 

*i(X) = (i.A)-MQf-QB(X)-QB(X,*o(X)) - 
-QB(*o (X),X)]. 

So, our function Ui{t) is related to an induced trajectory and, since this one 
is related to A^i, it is also related to this a.i.m. 

3. In the course of the numerical implementation of the method, qi will be 
computed in the points of the grid on [0, T], since it's values in these points will 
be used at the next level. 

5.3 Inductive definition of the high-order approximations 

Let us consider a k E N. We assume that for every < j < fc + 1, we already 
constructed Pj and q^ . Now, we define Pfe+2 as the solution of the problem 

Pfe+2-^Ap,+2+PB(pfe+2+q,+ i) = Pf, (30) 

P/c+2(0) = Puo, 

with qfc+i defined at the preceding step, and then set qfc+2 as 

qfe+2 = (i^A)"^ [Qf - QB(p^,^2) - QB(Pfc+2'qfe+i)~ 

-QB(qfc+i,Pfc+2) - QB(qfe,qfe) - q^] . (31) 

Naturally, the corresponding approximate solution of ([5])-® is defined by 

Ufe+2 (t) = Pk+2 {t) + qfc+2 (t) ■ 

Remarks. 

1. The right hand side of ((?T|l . that we denote by ^k+2{Pk+2 (t) , Qfc+i (t) , qfc (t)), 
is inspired from inductive the definition of the function qfe+2, m (t) of [15] . that 
is 

qfe+2,™ (t) = (!.A)"MQf-QB(p(t))- 
-QB (p(i) 7 qfe+l,m 

(t)) - QB (qfe+1,™ {t),p{t))~ 
-QB(qfe,™ {t))-ql„,{t)] . 
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Our functions Uk+2, k > are, in fact, approximations of the functions Ufc+2,Tn = 
p + qk+2.m that generate the induced trajectories in [15]. Our construction by- 
passes the construction of a.i.m.s. and is based directly upon that of the induced 
trajectories. We can call the sets {uk+2{t); t > 0} - approximate induced tra- 
jectories. 

2. The construction of the high accuracy a.i.m., Mk+2, is based upon 
the definition of the function qk+2,m of [15]. Aik+2 is the graph of ^k+2 ■ 
PTi -^Qn, 



*fc+2(X) = (i.A)-MQf-QB(X)- 

-QB (X, ^k+i (X)) - QB (X) , X) 

-QB(*fe (X))-D*fe (X)rfc(X)] 



(32) 



where D$fc (X) is the differential of €>fc (X), and 

Tk (X) = I. AX - PB (X + *fc (X)) + Pf. 

These a.i.m.s or some variant of these are used in the nonlinear Galerkin 
methods, and in the postprocessed high-order nonlinear Galerkin methods. 

3. If fc -|- 2 is the last level we construct, than we may compute cik+2 only at 
the moment of interest (T for example) as in the postprocessed method of [1]. 



6 Estimates of the error of the approximate so- 
lutions 

In the proof of the main result of this section, we need the following result that 
is a direct consequence of Lemma 1 from [4] . We denote by Vj j the coordinate 
of the function v with respect to the eigenfunction w* ; 

Lemma Let G(s) = ( i (s) ; ) and suppose that 



< cj for 0<j,l<m, l<i<A. 



The 



< 



E E 

j,k<.m 2—1 



(33) 



Now we can state and prove our main result. 

Theorem 2 The functions Uk (t) , k > 0, defined in the previous section, 
represent approximate solutions of the problem ([TJj-Q, and their accuracy in- 
creases with k. More precisely, the inequality : 
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|(u-u,)(t)| <C(5^/4+'=/2, 
holds for every k > and for t > {R) . 



(34) 



Proof We will prove our assertion by induction. 

1. We start with fc = 0. In [4^ the following estimate is proved, for 

\{p-Po){t)\<j^=C'S'/\ (35) 

where C is a constant, large for ly small. Actually, as can be seen from [1] 
this C is of the order of the product CqCi, with Cq, Ci the constants of our 

I Of I 

Theorem 1. Hence we can assume that C is of the form -^ ^23^3/ 2 for n great 
enough, with K a number depending on T but not on the data of the problem 
(see the proof of our Theorem 1). 

Let us observe that |po {t)\ is bounded for large times. Indeed, 

Ipo (01 = \p{t)+Po{t)-p{t)\<\p{t)\ + \po{t)-p{t)\ 
< Pq + C'S^^^ = i]q, for t large enough. 

The same observation is true for ||po {t)\\ and for |Apo {t)\ 

|Apo {t)\ < P2+ C'S^^^ ~ 772, for t large enough. 
In order to estimate the various norms of (q — qo) (t) we write 



q = (i^A)-i 



Qf-QB (p + q,p + q) + ^ 

at 



subtract from this relation the definition relation of qo, 
the obtained equality. In norm, we have 



(36) 

and apply i^A to 



|j.A(q-qo)| = QB(p + q)-QB(po) + ^ 

< |QB(p-po,p)| + |QB(po,p-Po)| + 
+ |QB(p, q)| + |QB(q, p)| + IQB (q, q) 

For the first term in the right side, with ([5]) we have: 



dq 



(37) 



|QB(p-Po,p)| < IIp-PoII IIpII 
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where we used once more the inequality L^S^^'^ < 1. Here and in the sequel, C 
denotes a generic constant (not depending on m but depending on i^, f , A). 

The same estimate holds for the second term. With ([7]), the third term 
yields: 



|QB(p,q)| < cilpp |Ap|^||q|| 

and the fourth 

|QB(q,p)| < cilqHAqpllpIl 

By using (|^D|) and all the above inequalities in ([57]) we obtain 

|^.A(q(i)-qo(i))| <C,5l/^ (38) 

for t great enough. As consequences 

||q(t)-q„(t)|| <C<5, |q(t)-qo(i)l<C5'/'. (39) 

Inequality (I35p and the second inequality above imply 

\u{t)-uo{t)\<C6'/\ (40) 

We must remark that, as is proved for u in |13j, we can prove that po is 
analytic in time, and more than that, it is the restriction of an analytic function 
of a complex variable to the real axis. This properties are transferred to q by 
its definition. Then, by using the Cauchy formula, it can be proved that 

|q'(t)-q;,(<)|<C<5^/^ (41) 



We also remark, for later use, that (|25p and the dissipativity of po imply 

|Aqo| < C, llqoll < CS^ and |qo| < CS for t > (R) ■ 

2. We now estimate |p — Pi| and |q — qi| . We have, by subtracting ((7f)) 
from (HSl) 



^(P-Pi) = i^A(p-pi)-PB(p + q-(pi+qo),p + q)- 

-PB (pi+qo, P + q (Pi+Qo)) 

From here, by using the semigroup of linear operators of infinitesimal generator 
uA, we obtain 



-PB (pi + qo, q - qo) - PB (pi + qo, p - pj} , 



:e''*^(p-Pi)(i) = e^*^{-PB(p-pi,u)-PB(q-qo,u)- 
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and, by integrating 



{p-Pi)W = 



(P-Pi)(0)- 



' f e-'^f*-^)^ {PB (p - pi, u) + PB (pi + qo, p - Pi)} ds - 
Jo 

f e-^(*-^)^ {PB (pi + qo, q - qo) + PB (q - qo, u)} ds. 
Jo 



Following [31 we use the inequalities pQ 
|A-''B(u,v)| < 
valid for 6 e (1/2,1) and [7] 



C|Ai-*u| |v| < C|Ai/2u| |v| 
C|u| lAi-'^vl < C|u| |Ai/2v| 



and obtain 



|(p - pi) {t)\ < |e--*A (p _ p^) (o)|+/J C{t~ sy' e-^(*-^) |(p - pi) {s)\ds+ 
+ |/o e-"-'*-^)^ [PB (pi + qo, q - qo) + PB (q - qo, p + q)] (s) ds . 

A form of Gronwall inequality ([7j, Lemma 7.1.1) implies 

l(p-Pi)WI< 

< C max^ |/o* e-'^(*-^)^ {PB (pi + qo, q - qo) + PB (q - qo, p + q)} (s) ds 

We must remark that the constant C above is of the order of . 
By using the method of [3] , we find the estimates for the coordinates of the 
several terms in the accolade: 



B(qo,q-qo),-fc < |qo|A2(q-qo) 



< CSS ^ cs^ 



B(q-qo,q)j-fe < |q-qol 
B(pi,q-qo)j-A; 



A2 q 



(42) 
(43) 



< 

< CS 



A5 (q - qo) (1(1 - Pm-,) P| + 1(1 - Pm-k) p|) 

1 1 



(44) 



B(q'^o,P),,fc < CKS'^'(-J- 



A5 

+ 1 m-fc+1 . 



where Pm-j represents the projection operator on the space spanned by the 
eigenfunctions corresponding to the eigenvalues in T„i^j and Xj = Xj^o. 
By using the inequalities ([331), (|42l), ((431) and 
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it follows that 



/* e-''(*-^)A {PB (qo, q - Qo) + PB (q - qg, q)} ds 
JO 



< CKS^ 



In order to estimate the term 
and the inequality 

j.k<m ^j,k m-j + l 



we use (IM 



< 



C 



(m + ir 
proved in ^3]. It follows 

Jo 



C5 



3/2 



The same estimate can be proved for 
hence finally we have 

|p-Pi|<C<S^/* 



jJg-.(t-.)A|pB (q_q„,p)}ds 

(45) 



We easily see that |pi| < yyg, ||pi|| < 77^, |Api| < 7^2. 

Now, in order to estimate the various norms of q — qi , we subtract (|28p 
from p6p . we apply the operator i/A, and take the norm in 7i of the resulted 
equality. After grouping the terms in a convenient way, we get 



|i/A(q-qi)| < |QB(p-pi,p)| + |QB(pi,p-pi)H 
+ IQB (p - Pi, q)| + IQB (qo, p - Pi) 
+ IQB (pi, q - qo)| + IQB (q - qo, p) 
dq 



IQB (q,q) 



dt 



(46) 



As we did for |q — qpl, we estimate one by one the terms from the right side. 
For the first one, we use ([8]) and the inequality < 1 : 

|QB(p-pi,p,)| < C2Ll/2||p_pi||||p|| 

< CL'/H^/*p^ < CS. 

The same estimate is valid for the second term. The third term is smaller 
than the first and for the fourth the following holds 



|QB(qo,p-Pi)| < cilqop lAqop IIp-PiI 
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For the two following terms we use ([5]) respectively ([TT|): 

|B(pi,q-qo)| < Cilpip |Apj5 ||q_qo|, 



|B(q-qo,p)| < C4|q-qo|^ ||q-qo||^ IIpII^ |Ap|^ 

The fifth and sixth terms are smaller than the first, respectively the second 
term, while for the seventh we have, with (fTTj) 

|B(q,q)| < C4|qN|qp||qp|Aq|5 
< G8. 

By using the above inequalities in (|46|) we obtain 

|j.A (q-qi)l < CL 

From here 

||q-qi|| <C5'/^ |q-qi|<CT2. (47) 

The arguments used to state the analyticity in time of qo remain valid for 
qi and the following relation follows 

This will be used later. By using (|45p and (|47| we now obtain 

|u-ui| < CW^/^ (48) 

3. The induction step. We assume that, for every < j < fc + 1 the 

inequalities 



p - 




< 




q- 




< 




|q' 




< 


^//j3/2+j72 






< 





hold. We prove that the above inequalities hold also for j = fc + 2 : 

(p - p,+2) (i) = e-''*^ (p-p,+2) (0)- 
- /o e-''^*-^)^ {PB (p - p,+2, u) + PB (pfc+2 + q^+i , p - p,+2) } ds- 
- /o e-''(*^^)^ {PB (pfc+2 + qfc+i , q - q^+i) + PB (q - q^+i , u) } ds. 



18 



As we did for \{p — p^) {t) \ , we obtain 

I (P - Pfc+2) {t)\ < |e-^*^ (p - Pfe+2) (0)|+/o Cit- s)-' e-^(*"^) |p - p,+,\ds+ 
+ |/o e--^^*-^)^ {PB {pk+2 + qfe+i, q - Qfe+i) + PB (q - q^^^, p + q)} ds 
The already cited Gronwall-type Lemma of [7] implies 
|(P - Pfe+2) it)\ < 

- ^o<t<T k" e"'^*"'^^ (pfc+2 + qfc+1, q - q,+i) + PB (q - q^^^, p + q)} 

We evaluate the coordinates of each term in the brackets following ^-'^i^-'^)-^: 



S(q-qfc+i,q] 



< |q-qfe+i 



Aaq 



S (Pn+2,q- q„+i)^.; < A2(q-q„+i) (|(I - P„_j) p„+2| + |(I - Pm-/) P«+2| 
The same arguments used for the terms involved in |(p — Pi) {t)\ lead to 



r e-''^*-^'^ {PB (q,+i, q - q,^,) + PB (q - q,+i, q) } ds 
Jo 

f\-.it-s)Ap^ (Pfe+2,q-qfe+i)d.s 
Jo 



= C'5('^+2)/2+5/4^ 



Analogously we can show that 



Jo 



< C(5(*=+2)/2+5/4^ 



and by putting these results together, it follows 

\ip~P,+,){t)\<CS^'^'^/'+'/\ (49) 
that confirms our induction hypothesis in what concerns p^ . It also follows that 

|Pfe+2|<%. l|Pfc+2|| < |Apfe+2|<?72- 

Now for |i/A (q — q_k^2) have 

i.A(q-qfc+2) = QB(p)-QB(p^^2) + QB(p,q)-QB(pfc^2,qfe+i) + 

+ QB(q, p)-QB(qfe^i, Pfe+2) + QB(q, q)-QB(q^, q^) + 

+q' q'fc 
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and 



|i/A(q-qfc+2)| < |QB(p-p;,+2,P)+QB(Pfe+2'P-Pfe+2)| + 

+ |QB(p - Pfc+2, q/c+i)+QB(qk+i,P - Pk+2))| + 

+ |QB(q- qfc+i,p)+QB(Pk+2,q- qk+i)| + 
+ |QB(q - q„ q)+QB(q„ q - q,)| + |q' - q'J ■ 

We can see, by using the induction hypothesis, and L}l'^8^l^ < 1, that, 
for the first two terms, we have 

|QB(p-p,^2,P)| I < c2Li/2||p-p,^Jh, 

|QB(Pfc+2,p-Pfe+2)| J - ' ^'+'11 

< CLi/2r;i,S5/4+(fc+i)/2 

The following term is smaller than the first. The fourth term can be estimated 
as follows 

I I i I I - II II 

|QB(qk+i:P-Pk+2)| < ci |qfc+i|2 |Aqfe+i|^ ||p-pfe+2|| 

For the fifth term, we obtain 

|QB(q-qfe+i,p)| < C4|q-q;,+i|^ ||q-q;,+i||* IIpP |Ap|5 

< Cp}/'p^/'(53/4+('=+i)/453/4+fc/4 

and for the sixth 



|QB(Pfc 



fc+2 



q- qfe+i)| 



< 
< 
< 



Ci |Pfc+2h |^Pfe+2|' II C 



qfc+i| 



Then, by using we obtain 
IQB(q- qfe,q 



|Aqp 



|QB(q„q-q,)| 



< C4|q-qfc+i|" ||q-qfc+i||" ||q|| 

< C'(5^/4+(fe+l)/4j3/4+/c/4^1/4 ^ ^^l + (fe+2)/2^ 

< C4|qfen|qfcp||q-qfcP|A(q-qfe)P 
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By using also the induction hypothesis on |q'— and by comparing the 
magnitude orders of the various terms, we find, successively. 



z/A(q- 


qfc+2)l 


< 


(7^1/2+(fe+2)/2 


l|q- 


- qfe+2|| 


< 




|q- 


- qfe+2| 


< 


(7^3/2+(fe+2)/2 



and these inequahties confirm our induction hypothesis. 
From (|49p and the above estimates it follows 

|u-u,+2| <C5^/^+('=+2)/2.n (50) 

7 Comments on the method 

1. A major advantage of our method is that we can use very low dimensional 
projection spaces for the approximations of p, since the accuracy of the approx- 
imate solution may be increased by using several iteration levels of the method. 

For example, if we choose m = 6, after having passed through five levels of 
the method we obtain an approximate solution U4(t) = \n{t) -t-q4(t) that bears 
an error of the order of 10"'^ since (5^^^ = ~ 0.0000032. At each level 

jj < j < 4 we will have to solve a system of 4x 36 + 4x6 = 168 ODEs, on 
the interval [0,T] in order to find the coordinates of Pj(i)i ^'^d to compute the 
coordinates of qj(t) by using algebraic relations. 

Here, as in all nonlinear Galerkin methods, problems appear due to f . If 
this function has a infinity of nonzero coefficients in its Fourier function, it will 
generate a infinite number of non-zero coordinates in q^ (t) . A truncation crite- 
rion must be applied and it will depend on f . Thus the number of coordinates 
of (\j [t) to be computed depends on j and on the given function f . 

If we chose m = 10, we need only four levels of the method for an error of 
the order of 10"^ (in this case <5"/^ = ~ 0.00000187). But at each level a 

number of 4 x 100 + 4 x 10 = 440 ODEs must be solved for the coordinates of 
Pj(i). Besides these, at each level j the coordinates of (\j{t) must be computed 
by algebraic relations resulted from the definitions. 

2. The program for the integration the systems of ODEs for Pj(t) should 
have the same structure for all j, only the coordinates of qj_i(t) remaining to 
be replaced in the nonlinear term. 

3. A comparison with the nonlinear Galerkin methods that use high- 
accurate a.i.m.s is necessary. 

The nonlinear Galerkin method based on the use of high accuracy a.i.m.s 
[5], [TT], applied to the Navier-Stokes problem and corresponding to our level 
k -\-2, k > 0, consists in solving the finite dimensional problem 

^-i.Ap+PB(p + *fc+i(p)) = Pf, (51) 

p(0) - Pu(0), 
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for the approximation jS of p = Pu. Here, as above, ^k+i '■ P'H Q'H, is the 
function defining an a.i.m. of high accuracy. The advantage of this method 
towards ours is that the system of equations for p is integrated only once. But 
the problem with solving (|5ip is that the definition of ^k+i requires those of all 
$j with j < k+1 and is very laborious (see Programming this must be very 
difficult. The structure of our method, with iterative levels, makes the compu- 
tations easier to program, and each level represents a certain approximation of 
the solution, so we can enjoy partial results. 

On another hand, all the computations for qo {t) , qi {t) , ...,q_k+i {t) (to- 
gether) in our method seem, at first glance, of the same order of complexity as 
those necessary for the evaluation of ^k+i (p (t)) in the course of the numerical 
integration of (|5ip in [2], [llj . However, in our method, a major simplifica- 
tion of the computations appears since ci'k_2 (from the definition of q^) may 
be approximated by the numerical derivative {qk-2(t) — qfc_2(i — h)) /h (since 
we have already computed qfc_2(t) at every time step). This must be com- 
pared with the definitions of z'j ,^ in or gj in [TT], that yield difficulties in 
the numerical integration programming. We must also remark here that the 
term D*fc-i (X) Tk-i (X) in the definition of *fc+i (X) (that is §^ with k + 1 
instead of fc -I- 2) is meant to approximate q'j,_]^ „j from the definition of q^+i^m- 
Hence, conceiving a method that uses directly the functions qfe+i,m instead of 
the a.i.m. s that are defined with the help of these functions is very natural. 

4. The memory of the computer is better organized in our method, since 
at the beginning of the computations for the level j we may erase from the 
memory the value of Pj-i{t) and keep only those of qj_i(t), qj_2(i)- 

5. In order to have not too many computations we may postprocess our 
solution pk+2{t) (at the last level) only at the end of the time interval [0,T], by 
adding qfc+2(r) = $fe+2 (pfc+2(r),qfe+i(T),qfc(r)), as is done in [U for ([51]). 

Acknowledgements 

Work supported by the Romanian Ministry of Education and Research via 
the grant CEEX-05-DII-25/5.11.2005 

References 

[1] p. Constantin, C. Foia§, Navier- Stokes Equations, Chicago Lectures in 
Math., Univ. of Chicago Press, IL, 1988; 

[2] C. Devulder, Martine Marion, A class of numerical algorithms for large 



time integration: the nonlinear Galerkin methods, SIAM J. Numer. Anal., 
29(1992), 462-483; 

[3] C. Foia§, O.Manley, R.Temam, Modelling of the interactions of the small 
and large eddies in two dimensional turbulent flows, Math. Modelling and 
Num. Anal, 22(1988), 93-114, ; 



22 



[4] B. Garcia- Archilla, Julia Novo, E.S. Titi, Postprocessing the Galerkin 
method: a novel approach to approximate inertial manifolds, SIAM J. 
Numer. Anal, 35(1998), 941-972; 

[5] B. Garcia-Archilla, Julia Novo, E.S. Titi, An approximate inertial mani- 
folds approach to postprocessing the Galerkin method for the Navier-Stokes 
equation, Mathematics of Computation, 68(1999), 893-911; 

[6] Adclina Gcorgescu, Hydrodynamic stability theory, Martinus Nijhoff Publ., 

Dordrecht, 1985; 

[7] D. Henry, Geometric theory of semilinear parabolic equations, Springer, 
Berlin, 1991; 

[8] M.S. Jolly, R. Rosa, R. Temam, Accurate computations on inertial mani- 
folds, SIAM J. Sci. Comp., 22(2001), 2216-2238; 

[9] Martinc Marion, R. Temam, Nonlinear Galerkin methods, SIAM J. Numer. 
Anal., 26(1989), 1139-1157; 

[10] Martine Marion, R. Temam, Nonlinear Galerkin methods: the finite ele- 
ment case, Numer. Math., 57(1990), 1-22; 

[11] Julia Novo, E.S. Titi, S. Wynne, Efficient methods using high accuracy 
approximate inertial manifolds, Numer. Math., 87(2001), 523-554; 

[12] J. C.Robinson, Infinite- dimensional dynamical systems; An introduction to 
dissipative parabolic PDEs and the theory of global attractors, Cambridge 
University Press, 2001; 

[13] R. Temam, Namier- Stokes equations and nonlinear functional analysis, 
CBMS-NSF Reg. Conf. Ser. in Appl. Math., SIAM, 1995; 

[14] R. Temam, Infinite- dimensional dynamical systems in mechanics and 
physics. Springer Appl. Math. Sci., 68, Springer, New- York, 1997; 

[15] R. Temam, Induced trajectories and approximate inertial manifolds. Math. 
Mod. Num. Anal., 23(1989), 541-561; 

[16] R. Temam, Attractors for the Navier-Stokes equations, localization and 
approximation, J. Fac. Sci. Univ. Tokyo, Soc. lA, Math., 36(1989), 629- 
647. 



23 



